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ABSTRACT 

The hyperluminous X-ray source HLX-1 in the galaxy ESO 243-49, currently the best intermediate-mass black 
hole (BH) candidate, displays spectral transitions similar to those observed in Galactic BH binaries, but with a 
luminosity 100-1000 times higher. We investigated the X-ray properties of this unique source by fitting multi- 
epoch data collected by Swift, XMM-Newton, and Chandra with a disk model computing spectra for a wide range 
of sub- and super-Eddington accretion rates assuming a non-spinning BH and a face-on disk (i = 0°). Under these 
assumptions we find that the BH in HLX-1 is in the intermediate-mass range (~2 x 10 4 M 0 ) and the accretion 
flow is in the sub-Eddington regime. The disk radiation efficiency is rj = 0.11 ± 0.03. We also show that the 
source does follow the Lx oc 7’ 4 relation for our mass estimate. At the outburst peaks, the source radiates near 
the Eddington limit. The accretion rate then stays constant around 4 x 10~ 4 M 0 yr _1 for several days and then 
decreases exponentially. Such “plateaus” in the accretion rate could be evidence that enhanced mass-transfer rate 
is the driving outburst mechanism in HLX-1. We also report on the new outburst observed in 2011 August by the 
Swift X-Ray Telescope. The time of this new outburst further strengthens the ~1 year recurrence timescale. 

Key words: accretion, accretion disks - black hole physics - galaxies: individual (ESO 243-49) - methods: data 
analysis - X-rays: individuals (HLX-1) 
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1. INTRODUCTION 

Ultraluminous X-ray sources (ULXs) are defined as off- 

nucleus extragalactic sources showing X-ray luminosity exceed- 

ing 3 x 10 39 erg s' 1 assuming isotropic emission (see Roberts 

2007). Even if their nature is still in dispute, it is likely that 

their huge luminosity is produced by accretion of matter onto 

a black hole (BH). Three explanations for their nature have 
been considered. (1) ULXs may be X-ray stellar mass BH 
binaries (BHBs) similar to those observed in our Galaxy, but in 

a more extreme version of the very high state (e.g., Remillard & 
McClintock 2006), the ultraluminous or wind-dominated state 
(e.g., Gladstone et al. 2009). In such states, the source would 
be able to radiate above the Eddington limit. It is still to 
be understood why this ultraluminous state is so rarely ob- 

served in X-ray binaries. (2) King et al. (2001) proposed that 
the emission of ULXs is highly anisotropic. In this case, the 
requirement to have super-Eddington emission is alleviated. 
The nature of this anisotropy could be due to either geo- 
metrically thick disks funneling the X-ray photons produced 
in the inner parts of the accretion disks (King 2009) or due 
to relativistic beaming of a jet or strong outflow. However, 
the discovery of several optical and radio nebulae around 
ULXs (e.g., Pakull & Grise 2008) and quasi-periodic oscil- 
lations in M82 X-l (Strohmayer & Mushotzky 2003) argue 
against strong beaming. (3) The final and the most exciting 
explanation is that some ULXs are accreting intermediate- 
mass BHs (IMBHs) with masses ranging from ~1OOM 0 to 


~10 5 Mq (e.g., Colbert & Mushotzky 1999). The existence 
of such IMBHs will naturally alleviate the need for super- 
Eddington emission. It is clear now that the IMBH interpreta- 
tion is not valid to explain the ULX population as a whole (see, 
e.g., Roberts 2007). However, the most luminous ULXs, the so- 
called hyperluminous X-ray sources (HLXs; Gao et al. 2003) 
with X-ray luminosities above 10 41 erg s _1 , are good candi- 
dates. Finding convincing evidence for the existence of IMBHs 
is important for the growth of supermassive BHs via mergers or 
accretion episodes (e.g., Micic et al. 2007), dark matter studies 
(e.g., Fornasa & Bertone 2008, but see also Bringmann et al. 
2009), cosmology (e.g.. Trend & Stiavelli 2007), and gravita- 
tional wave detection (e.g., Matsubayashi et al. 2004; Amaro- 
Seoane & Santamarfa 2010). 

Farrell et al. (2009) reported the serendipitous discovery of a 
ULX candidate 2XMM JO 1 1 028 . 1 -46042 1 , referred to hereafter 
as HLX- 1 , located in the outskirts of the edge-on spiral galaxy 
ESO 243-49 at a redshift of 0.0224 (Wiersema et al. 2010). 
From its maximum luminosity reaching ~1.3 x 10 42 erg s _1 
at peak and assuming that the source luminosity reached up 
to 10 times the Eddington limit (indd), Farrell et al. (2009) 
derived a BH mass of more than 500 M 0 . HLX-1 is so far 
the best candidate to harbor an IMBH. Apart from its extreme 
luminosity, which has been observed many times over the past 
3 years by different X-ray satellites (Swift, XMM-Newton, and 
Chandra), HLX-1 is unique among other ULXs because it 
is the only one for which clear spectral hysteresis similar to 
those observed in Galactic BHBs (GBHBs) are seen. Indeed, 
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Figure 1. Swift-XKY photon counting 0.3-10 keV light curve of HLX-1 up to 2011 September 5. The Swift-XRT light curve was obtained using the Swift — XRT 
light-curve generator web interface with a binning of at least 60 counts per bin (Evans et al. 2009). Left: the vertical dotted lines mark the times of the XMM -Newton 
and Chandra observations. Right: close-up view of the X-ray light curve around the outburst peaks PI, P2, and P3. The zero time corresponds to the time when we 
first observed a re-brightening from the source (PI = 2009 August 16 and P2 = 2010 August 29) and the peak of the outburst P3 (201 1 August 15). The square on 
the middle panel is the predicted on-axis Chandra count rate into the 0.3-10 keV Swift-XRT energy band. Bottom panel: the thin points were obtained using a lower 
temporal binning for data from MJD = 55,788 to MJD = 55,790. 

(A color version of this figure is available in the online journal.) 


Servillat et al. (201 1), using XMM-Newton, Chandra, and Swift 
data, confirmed the results presented in Godet et al. (2009) 
that HLX-1 underwent transitions from the high/soft state to 
the low/hard state. The Swift X-Ray Telescope (XRT) light 
curve so far covers four outbursts from 2008 to 201 1 (hereafter 
P0, PI, P2, and P3 — see Figure 1). While we had only one 
observation in 2008 October during the P0 outburst, the PI and 
P2 outbursts display well-sampled fast-rise/exponential-decay- 
like temporal profiles. Thanks to our dedicated Swift-XRT target 
of opportunity, we caught the rise of the P3 outburst from 2011 
August 15 (Godet et al. 2011). A precursor could be seen prior 
to the P3 outburst peak in Figure 1 . Two re-flare events are also 
visible during the PI and P3 outbursts. The outbursts appear to 
be separated by a recurrence timescale of nearly a year. Recently, 
Lasota et al. (2011) interpreted the X-ray light curve as the result 
of enhanced mass-transfer rate onto a pre-existing accretion disk 
around an 1MBH, when an asymptotic giant branch star orbiting 
along an eccentric orbit with a period of ~1 year is tidally 
stripped near periastron. 

In order to further investigate the nature of this unique source, 
it is essential to put some constraints on the accretion flow and 
the BH mass. Lower and upper limits from radio, optical, and 
X-ray observations have been derived (see Wiersema et al. 2010; 
Servillat et al. 2011; Webb et al. 2012). Davis et al. (2011), 
using their relativistic accretion disk model bhspec to fit X-ray 
spectra when the source was in various spectral states, put some 
constraints on the BH mass. They showed that the inclination (;) 
has a strong influence on both the mass and the spin (a*) of the 
BH derived by the model. Due to degeneracies in their best-fit 
parameters, they were only able to give a range of possible BH 
masses from 3000 M 0 (where the limit corresponds to i = 0°, 
(L x /L E a d) = 0.7, and a* = — 1) to 3 x 10 3 M 0 (where the 
limit corresponds to i = 90° and a* = 0.99). However, the 
inclination is likely to be less than 60° -70° due to the lack of 
observed eclipse dips assuming HLX-1 is a binary system. These 
results were obtained neglecting the effects of radial advection 
that are important at luminosities above the Eddington limit L E dd 
with L Ed d = L3 x 10 38 (M/M q ) erg s -1 . Radial advection is a 


mechanism enabled to stabilize disks in the super-Eddington 
accretion regime (see, e.g., Abramowicz et al. 1988). This 
regime is reached when the accretion rate (M) is larger than 
M E dd = i 1 ~ l (L E dd/c 2 ) with c being the speed of light and ii 
being the radiation efficiency. In this regime, the vertical disk 
structure could be geometrically thick (H/R > 0.5 with Ii and R 
being the scale height and the radius of the disk, respectively). 
However, for a range of moderate super-Eddington M-values 
there are still disk solutions for which the vertical disk structure 
can be considered as relatively thin. They are often referred 
to as slim disk solutions (Abramowicz et al. 1988). For large 
accretion rates (m ^ 10 with m = {Me 2 / L E dd)), electron 
scattering (opacity and Comptonization) has several effects on 
the emergent spectra: (1) they can be highly distorted and no 
longer look like a multi-color blackbody (BB) spectrum (e.g., 
Kawaguchi 2003, hereafter K03); (2) the LocT 4 relation with 
T being the color temperature at the inner disk radius for a 
standard ex - cl i s k (Shakura & Sunyaev 1973) is no longer valid 
(see Fukue 2000; K03). In addition, slim disks may extend 
beyond the innermost stable circular orbit (ISCO) due to non- 
negligible pressure support resulting in higher disk temperatures 
(e.g., Abramowicz et al. 1988, 2010; Watarai et al. 2000). 

In this paper, we make use of the K03 disk model in order to: 
(1) investigate the accretion flow and BH properties in HLX-1 
through spectral fitting of a multi-epoch and multi-instrument 
data set using a model that computes a wide range of sub- 
and super-Eddington accretion disk solutions; (2) investigate 
whether radial advection and electron scattering have important 
effects on the emergent spectra; (3) compare our results with 
those performed on less luminous ULXs using the same model 
that always favored critical or super-Eddington accretion onto a 
stellar mass BH (e.g., Yoshida et al. 2010; Foschini et al. 2006; 
Vierdayanti et al. 2006; Okajima et al. 2006). 

The paper is organized as follows. In Section 2, we describe 
the data used to perform the spectral analysis together with the 
data reduction. Section 3 presents some possible observational 
evidence for effects of radial advection and/or electron scatter- 
ing as well as the detection of possible X-ray lines. In Section 4, 
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Table 1 

Log of the XMM-Newton, Chandra, and Swift Observations 


Obs. Name 

Instrument 

Obs. ID 

Start Date 

End Date 

Good Time 
(ks) 

XMMl 

XMM-Newton EPIC-pn 

0204540201 

2004 Nov 23 

2004 Nov 23 

22 

XMM2 

XMM-Newton EPIC-pn, 

0560180901 

2008 Nov 28 

2008 Nov 28 

51 

XMM2 

XMM-Newton RGS1, RGS2 

0560180901 

2008 Nov 28 

2008 Nov 28 

51 

XMM3 

XMM-Newton EPIC-pn 

0655510201 

2010 May 14 

2010 May 15 

62 

Chandra 

Chandra ACIS 

13122 

2010 Sep 6 

2010 Sep 7 

10 


Swift-XRT 

00031287(001-139) 

2008 Oct 24 

2011 Oct 10 

360 


we present the fitting results obtained with the K03 disk model. 
In Section 5, we discuss the implications of the fitting results 
on the BH mass, the accretion rate, and the disk structure. The 
main conclusions are given in Section 6. 

2. OBSERVATION LOG AND DATA REDUCTION 

We consider X-ray data coming from three observatories: 
XMM-Newton, Chandra (ACIS), and Swift. The data used in 
this study cover the different spectral states observed dur- 
ing the outbursts of HLX-1 (see Table 1). We used the same 
XMM-Newton and Chandra spectra as those presented in 
Servillat et al. (201 1). The nomenclature used here is the same 
as that in Servillat et al. (2011) and Farrell et al. (2009) for the 
XMM-Newton data (see Table 1). For the XMM-Newton data, we 
only used the EPIC-pn data because they offer the best statistics 
for a given observation. Please refer to Servillat et al. (2011) for 
the details of the data reduction on the XMM-Newton and Chan- 
dra data. In order to investigate for the presence of X-ray lines 
during the outbursts, we also used the Reflection Grating Spec- 
trometer (RGS) data from the XMM2 observation (see Table 1). 
The RGS data have been reduced using the standard XMM- 
Newton Science Analysis Software task rgsproc and the most 
recent calibration files released in 2011 February. For the spec- 
tral analysis, we only considered the first order of the RGS data. 

All the Swift-X RT photon counting data were processed us- 
ing HEASOFT v6. 1 1 and the latest calibration files (CALDB 
version 3.8). This new CALDB includes a new gain file enabled 
to correct the data for charge traps that accumulate on the CCD 
due to radiation damage (see Pagani et al. 2011). These traps 
can induce some energy offsets. The data were processed us- 
ing the tool xrtpipeline vO.12.6. We inspected all segments 
of data to search for epochs when the background was en- 
hanced by the presence of hot pixels induced by a high CCD 
temperature ( T > —55°). For these time intervals, we checked 
whether there were some hot pixels in the extraction regions 
used to extract the spectra. Thus, we excluded the data from 
segments 00031287125 (0.8 ks) and 00031287129 (0.2 ks) 
because the data were severely contaminated by hot pixels due 
to a high ( T > — 51°C) CCD temperature. We also excluded 
the first orbit (~0.6 ks) of segment 00031287023 for the same 
reason. These latter bad data induce an artificially high count 
rate in the X-ray light curve at MJD = 55,209 (during the 
PI outburst — see Figure 1). However, the re-flare event seen 
around that time is not due to hot pixel contamination. We used 
the grade 0-12 events, giving slightly higher effective area at 
higher energies than the grade 0 events, and a 20 pixel (47.2 arc- 
sec) radius circle to extract the source and background spectra 
using xselect v2.4b. The background extraction region was 
chosen in order to be close to the source extraction region 
and not to contain any XMM-Newton sources. The ancillary 


response files were generated using xrtmkarf v0.5.9 and expo- 
sure maps generated by xrtexpomap v0.2.7. The response file 
swxpcOto 1 2s6_ 20010101v01 2 .rmf is used to fit the spectra. To 
extract the Swift-X RT spectra, we divided the data over different 
intensity ranges for the outbursts P0, PI, and P2. As the work on 
the first three outbursts had been completed before the most re- 
cent outburst P3, we decided to compare the latest outburst with 
the previous three by dividing the data over time instead. Table 2 
summarizes the nomenclature used for the Swift-X RT data. 

3. SPECTRAL ANALYSIS 

All the spectra were grouped at a minimum of 20 counts per 
bin to provide sufficient statistics to use the yf minimization 
technique. The XMM3, Swift-X RT, and Chandra spectra were 
fitted in the 0.3-10 keV energy range within xspec v12.7.0 
(Arnaud 1996), while the XMM1 and XMM2 spectra were 
fitted in the 0.2-10 keV energy range. The Chandra data, 
being moderately piled up, were fitted using the pileup model 
(Davis 2001) within xspec with a frame time fixed at 0.8 s 
(see Servillat et al. 2011 for more details). For each model, 
we determined the appropriate value of the grade morphing 
parameter. For the data used in this paper, the statistical errors 
are dominant over the instrumental systematics. The RGS1 
and RGS2 spectra were binned at twice the resolution of the 
instrument (Al = 0.2 A). Since with this choice of binning 
there are fewer than 20 counts per resolution bin, we used 
the C-statistic (Cash 1979) available within xspec for the 
spectral fit. 

The total column density was fixed at the best-constrained 
value of A H = 4 x 10 20 cm -2 from the XMM2 observation 
(see Farrell et al. 2009). The Galactic absorption column in 
the direction of the source is equal to 1.8 x 10 20 cm -2 (e.g., 
Kalberla et al. 2005). For all fits, the absorption is modeled 
using the wabs photoelectric absorption model. With a source 
redshift of 0.0224 (Wiersema et al. 2010), we adopted a source 
distance of di = 95 Mpc using the cosmological parameters 
from the WMAP5 results (Hq = 71 km s 1 Mpc -1 , G M = 0.27, 
and £l\ = 0.73). All the errors quoted below are given at 
a 90% confidence level for one parameter of interest (i.e.. 
Ax 2 = 2.706). 

3.1. Hardness-Intensity Diagram 

We fitted the Swift-X RT spectrum using a multi-color disk 
diskbb and/or powerlaw model. This phenomenological 
model is often used to fit the spectra of GBHBs. In that context, 
the power-law component is often interpreted as emission 
produced by a corona of hot electrons located in the inner 
regions of the accretion disk, while the diskbb component 
is interpreted as the emission from the accretion disk 
(see Remillard & McClintock 2006). The photon index of the 
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Table 2 

Nomenclature of the Swryf-XRT Spectra 


Spectrum Name 

Comment 

Sip 

Time when the source reached its maximum luminosity 
for the first outburst PI (a count rate of CR > 0.02 counts s -1 ) 

Sip 

Same but for the second outburst P2 (CR > 0.02 counts s -1 ) 

S\p+2p 

Sum of the S\ p and S 2 P spectra 

^after2 p 

~7 days after the peak of the second outburst P2 (CR > 0.02 counts s -1 ) 

ScRO 

0.01 < CR < 0.02 counts s _1 

ScRl 

0.015 < CR < 0.02 counts s -1 

ScR2 

0.01 < CR < 0.015 counts s _1 

ScR3 

0.004 < CR < 0.01 counts s -1 

ScR4 

0.002 < CR < 0.004 counts s -1 

ScR5 

CR < 0.002 counts s^ 1 

P3~ 

Outburst in 201 1 Aug (55,788 < MJD < 55,809) 

P3+ 

Outburst in 201 1 Aug (55,710 < MJD < 55,844) 

S,i 

55,788 < MJD <55,802 

S, 2 

55,803 < MJD <55,822 

S, 3 

Re-flare event in the P3 outburst (55,826 < MJD < 55,844) 


Table 3 

Summary of the Spectral Parameters when Fitting the Swift-XRT Data Using the diskbb Model 


(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

GO) 

Name 3 

kT 

Nbb 

r 

Nr 

HR 

f^tot 

f-'disk 

log(Lboi) 

X 2 /dof 

Si p+2 p 

0.24 ±0.01 

20+ ® 



0.06 ±0.01 

12.0 ±0.4 

12.0 ±0.4 

42.2 ± 0.03 

42/36 

Sip 

0.24 ± 0.02 

20+6 



0.06 ±0.01 

12.0 ±0.5 

12.0 ±0.5 

47 i q+ 0.04 
^■ iy -0.05 

18/21 

Sip 

0.22 ± 0.02 

30/\ 4 0 



0 04 +0 - 01 
0.02 

13.3 ±0.7 

13.3 ±0.7 

47 77+0.05 

h-z.zj_oo 6 

24/15 

5"after2 p 

0.18 ±0.02 

71+58 
71 -33 



0 022 +0 006 
U.UZZ_Q QOS 

13.3 ±0.8 

13.3 ±0.8 

42.25 ± 0.07 

54/67 

ScRO 

n 1 O+0.02 

u - i8 -o.oi 

^6 +14 

JO-12 

9 9+0.5 
Z ’ Z -0.6 

9 Q+ 1 - 8 
z "-1.4 

0 2O +0 04 
u - zu -0.06 

H c+0.3 
'■ J -0.6 

5 6 +0 - 3 

41.90 ±0.04 

54/67 

ScR3 

0.17 ±0.05 

11 + 36 
1 —8 

9 9+0.5 
Z ‘ Z -0.6 

1 .8 ± 1.5 

0.37 ±0.15 

73 73 + 0.4 
J - 0.8 

2 0 +0 - 3 
z>u -0.7 

41.46+ 0 o 1 1 °3 

54/67 

ScR4 



9 9+0.5 
Z ’ Z -0.6 

2.0 ±0.7 

1 2 +0 - 9 
i- — 0.5 

1.3 ±0.2 



54/67 

ScR5 



9 9+0.5 
Z ’ Z -0.6 

0 3 +0 - 3 
u -0.2 

1.2 ±0.9 

0.20 ± 0.06 



54/67 

ScRl 

0.18 ±0.02 

41 -19 

2.2 b 

2.9 ± 1.2 

o 1 5+0-08 
u,iJ -0.06 

o 7+0.9 
*•'-1.6 

7.8+°i' 5 

42.03 ± 0.06 

45/51 

ScR2 

0.17 ±0.02 

40+ 2 t7 

2.2 b 

3.0+0.7 

o 99+O.O6 
-0.07 

6 7 +0 - 7 
°- / -0.8 

4 8 +0 - 3 

41.85 ±0.06 

45/51 

S,1 

0.20 ± 0.02 

73 0+25 
*"-15 



0 029 +0 012 

U.UZy_ooio 

12.4 ±2.0 

12.4 ±2.0 

42.22 ± 0.07 

15.2/15 

S ,2 

0.20 ± 0.02 

41-12 



0 023 +0 006 

G.uz,j_0008 

10.5 ± 1.0 

10.5 ± 1.0 

42.15 ±0.05 

17.2/22 

S, 3 

0.18 ±0.03 

44+25 

2 .2 b 

9 0+1.9 
Z -°-2.0 

pi i q+0.13 
U,18 -0.08 

8.6 ± 1.7 

6.8 ± 1.2 

41.98 ±0.10 

9.5/7 


Notes. Columns: (1) spectrum name; (2) inner disk temperature in units of keV; (3) the blackbody normalization; (4) photon index of the power-law component; 
(5) normalization of the power law in units of 10~ 5 photons kcV 1 cm~ 2 s -1 ; (6) hardness ratio defined as the ratio of the observed 1.5-10 keV flux over 
the observed 0.3-10 keV flux; (7) unabsorbed 0.2-10 keV total luminosity in units of 10 41 erg s 1 ; (8) unabsorbed 0.2—10 keV disk luminosity in units of 
10 41 erg s -1 ; (9) logarithm of the unabsorbed 0.01-20 keV disk luminosity in units of 10 41 erg s -1 ; (10) / 2 -value and number of degrees of freedom. 
a See Table 2 for the nomenclature of the different spectra. 
b Fixed parameter. 


power law was tied together for the Scro, i,2,3,4,5 spectra because 
we found F -values consistent within the errors when fitting 
the spectra individually. The power-law normalization and 
the parameters of the diskbb model were left free to vary 
independently between spectra. The S/1,2 spectra are well 
fitted by an absorbed diskbb model. We did not obtain a good fit 
(x 2 /dof = 14.5/8) using an absorbed diskbb model for the S , 3 
spectrum. There are residuals left at high energies. When adding 
a power-law component with T free, the T parameter was not 
constrained. So, we decided to freeze this parameter at the same 
value as found for the Sato, 1,2, 3, 4, 5 spectra (T = 2.2). The fit is 
improved (A/ 2 = 5 for 1 dof, i.e., an improvement at a >2 a 
significance level). For each spectrum, we then computed the 
hardness ratio (HR) of the 0.3-1. 5 keV observed flux over the 
1.5-10 ke V flux. This is different from the Swift HR presented in 


Godet et al. (2009) and Servillat et al. (201 1) that was based on 
count rate. The HR values, the best-fit spectral parameters, and 
the 0.2-10 keV unabsorbed luminosity ( L = 4 ltd \ F with F the 
0.2-10 keV unabsorbed flux) for each spectrum are summarized 
in Table 3. 

Figure 2 shows the Swift-X RT HR versus the 0.2-10 keV 
unabsorbed luminosity. The points from the Chandra ( kT = 
0.21 ± 0.01 keV and N = 37+“), XMM2 (kT = 0.17 ± 
0.01 keV, N = 40+ '|q, T = 2.1 ± 0.4, and N r = 1.9+° 0 8 6 x 1(T 5 
photons cm~ 2 s -1 keV” 1 ) and XMM3 10 (T = 2.2+gg and 


10 We also added a mekal component with kT = 0.44+ q 1 ^ keV to take a 
possible contamination from the galaxy into account (see Servillat et al. 2011). 
Only the power-law component was used to compute the hardness ratio. 
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Figure 2. Swift-X RT hardness vs. intensity diagram. The XMM-Newton 
(triangles) and Chandra (square) points are also reported on the diagram. The 
arrows indicate the evolution from the peak to the end of the outbursts. 

(A color version of this figure is available in the online journal.) 
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Nr = 2.9 ± 0.7 x 10~ 6 photons cm -2 s^ 1 keV” 1 ) observations 
are also reported on the plot. 

Figure 2 shows that we have a good agreement be- 
tween the different instruments. The HLX-1 track in the 
hardness-intensity diagram (HID) is clearly reminiscent of 
that observed in GBHBs. Between S i p + 2 p and .S' a ftcr 2 /), there 
is a softening of the source at nearly constant luminosity 
(~ 1 .2 x 10 42 erg s _1 ). Such a track in the HID is observed in GB- 
HBs (see, e.g., Meyer-Hofmeister et al. 2009). However, the sit- 
uation appears to be different because the disk temperature near 
the peak of the outbursts PI, P2, and P3 varies at nearly constant 
luminosity (from 0.24 ± 0.02 keV to 0.18 ± 0.02 keV — see 
Table 3). Note that if we consider the 3 a errors on the tem- 
perature (i.e.. A/ 2 = 9 for one parameter of interest), then 
we find that all data have consistent kT -values given the lim- 
ited statistical quality of the X-ray data used. However, when 
merging the S\ p and Sip spectra ( S\ P + 2 P ) and the S t \ and S, 2 
spectra (5 1 , 1+f2 ), the fcT-values (kT = 0.240 4 qq 23 0 for S\ p+ 2 P 
and kT = O^O 4 ^ for S, i+, 2 ) are no longer consistent within 
the 3 a errors. 


3.2. Investigating the L-T Relation 

Figure 3 displays the disk temperature versus the unabsorbed 
0.2-10 keV luminosity. We fitted our L-kT points using a 
La.T a relation, where a was left as a free parameter. We 
obtained a good fit with log(L) = (2.4 ± 0.7) log T + 43.7 ± 0.5 
(la errors). Feng & Kaaret (2007) showed that the luminosity 
L varies as L ocT~ xl±0 - 5 in the ULX NGC 1313 X-2. Using 
a larger sample of ULXs, Kajava & Poutanen (2008) found 
a La T~ a relation. King (2009) suggested that the L oc T 4 
relation is expected in the super-Eddington regime provided 
that the geometrical beaming factor varies as b(xm~ 2 . From 
Figure 3, the observed trend in the L-kT space does not follow 
such a correlation. Mizuno et al. (2001) found that some ULXs 
do present a L oc T 2 relation. They interpreted such a correlation 
as a possible signature for moderate super-Eddington accretion 
(see also Fukue 2000; K03). However, we cannot exclude a T 4 
relationship at the 3a level (see also Servillat et al. 2011). In 
Figure 3, we also show the best fit obtained using the L ocT 4 
relation. 



Figure 3. Evolution of the unabsorbed 0.01-20 keV disk luminosity with the 
disk temperature. The values for Swift (diamonds), XMM-Newton (XMM2: 
triangles in orange), and Chandra (square in blue) are shown on the plot. The 
thick dot-dashed line corresponds to the best fit obtained using the LocT 4 
relation, while the thick solid line corresponds to the best fit obtained using the 
L a T“ relation with a ~ 2.4. The thick dashed line corresponds to the L tx T 4 
prediction in the super-Eddington regime (Fukue 2000; King 2009). 

(A color version of this figure is available in the online journal.) 

3.3. Evidence for Possible Line Features 

When fitting the S\ P + 2 P spectrum, we found some strong 
residuals around 0.6 keV (see the top left panel in Figure 4). 
There may be some lower significance residuals at higher energy 
as well. We separately investigated the S\ p and St p spectra, and 
both show some residuals around 0.6 keV, in particular for S 2 P . 
Given that the residuals are close to the instrumental oxygen 
edge, we investigated if the residual may have an instrumental 
origin. We extracted a background spectrum using different parts 
of the field of view and encompassing the time when the source 
was at the peak. The fit of the background spectrum does not 
reveal any residuals around the O edge. We checked the images 
and found no evidence of contamination by hot pixels. When the 
observations were performed the source location on the CCD 
was not close to the bad columns. The application of the charge 
trap correction did not change the residuals. So, the residuals are 
unlikely to be instrumental. We added a Gaussian line redshifted 
at the galaxy distance ( z = 0.0224) with an intrinsic width equal 
to zero, because this parameter is otherwise not constrained. 
The fit was improved by A/ 2 = 12.7 for 2 dof. The line 
centroid is equal to 0.6 18 4 0042 keV and the equivalent width 
is EW = 93 47 5 5 , eV (see the top right panel in Figure 4). The 
disk parameters do not change with kT = 0.25 ± 0.01 keV and 
N = 1 6 + 4 (see Table 3 for a comparison). In the rest of the paper, 
we did not include the Gaussian line when fitting the S\ P+ 2 P 
spectrum because it did not change the spectral parameters of 
the K03 disk model considered in Section 4 significantly. 

In order to further investigate the soft X-ray emission at 
~0.6 keV, we inspected the RGS data available from the XMM2 
observation. Because the RGS1 and 2 spectra are background 
dominated below 0.5 ke V and above 1 .3 ke V, we only considered 
this energy range in our analysis. In the 0.5-1 .3 keV energy band 
a total of 1230 and 1170 counts were collected in the RGS1 
and RGS2, respectively (with 520 and 605 net counts). We first 
used an absorbed diskbb plus a power-law component to fit 
the data, fixing the photon index of the power-law component 
and the absorbing column density to the values derived from 
the XMM2 data (Farrell et al. 2009; Servillat et al. 2011). We 
found a temperature (kT = 0.19 ± 0.02 keV) consistent with 


5 


The Astrophysical Journal, 752:34 (12pp), 2012 June 10 


Godet et al. 



Rest Energy (keV) 

Figure 4. X-ray lines present in the Swift— X RT and XMM-Newton RGS data. Top: Swift— X RT spectrum when the source luminosity was at peak for the first and second 
outbursts (S\p+ 2 p). Left: using only an absorbed diskbb model. Right: same model with the addition of a Gaussian line redshifted at the galaxy redshift (z = 0.0224). 
The dotted lines correspond to the different components of the model. Bottom: ratio of the RGS 1 (red) and RGS2 (black) spectra from the XMM2 observation over 
the model (an absorbed diskbb + power-law model). The vertical dashed lines correspond to the centroids of the detected emission lines. 

(A color version of this figure is available in the online journal.) 


that derived from the EPIC data (see Table 3; see also Farrell 
et al. 2009; Servillat et al. 201 1). Although the fit is statistically 
acceptable (C-statistic = 160.4 for 113 bins) a first inspection 
of the RGS1 and 2 residuals (see the bottom panel in Figure 4) 
reveals that features are present at ~0.6 keV, ~0.9 keV, and 
~1 keV. 

We then added to the baseline continuum three unresolved 
Gaussian emission lines redshifted at the galaxy distance. The 
improvement in the fit due to the addition of the ~0.6 keV 
emission line is AC /dof = 26.2/2 (corresponding to a detection 
significance >99.9% confidence). The emission line is detected 
at an energy centroid of E = 0.642 ± 0.002 keV, a flux of F = 
(1.91 ± 0.6) x HI ' photons cm' 2 s~ 1 , and an equivalent width 
of EW = 20 ± 7 eV. The energy centroid is consistent with that 
derived from the Swift-XRT data, and the line equivalent width 
seems to decrease with the luminosity of the continuum. An 
additional line is detected at E = 0.880^0002 keV (AC/dof = 
12.2/2) with a flux of F = (0.93±0.6)x 10“ 5 photons cm' 2 s -1 
(EW = 27 ± 12 eV). Finally, the addition of the third Gaussian 
line at E = 0.98 ± O.OlkeV (EW = 31 ± 22 eV) gave a 
much lower statistical significance (AC/dof = 7.0/2). Even 
if additional data are needed to further investigate the origin of 
these lines, we speculate that the most likely identification for the 
~0.6 keV and ~0.9 keV lines is O vm Lyot and Fe xvin-Fe xix. 
The third line, if any, might be associated with Ne x Lya. 

4. INVESTIGATING THE ACCRETION 
DISK STRUCTURE 

4.1. A Simple Approach 

In order to investigate whether the Lx oc T~ 1A relation 
we found in Section 3.2 could be the result of a slim disk 
surrounding the BH, we first used the diskpbb model instead of 
the diskbb model. The diskpbb model is a multi-temperature BB 


Table 4 

Summary of the p - value Derived Using the diskpbb Model for the 
XMM-Newton , Swift-XRT, and Chandra Data 


Spectrum 

P 

XMMl 

0 < 0.54 

XMM2 

0 63 +0,37 

Chandra a 

n 57+0.33 
u -J/_ o.o7 

Sip 

0.71 ±0.2 

Sip 

0.53 - 1 

5"after2 p 

o 55+0.45 
U - J - 0.05 

ScRO 

0 7 +0 - 3 
u>/ -0.2 


Note. a The value was obtained fixing Nu to 4 x 
10 20 cm' 2 . 


disk model where local disk temperature is given by T (r) oc r~ p 
with p a free parameter (e.g., Mineshige et al. 1994; Hirano 
et al. 1995; Watarai et al. 2000). When p = 0.75, the model 
is equivalent to a diskbb model. For p < 0.75, radial advection 
starts to become important (e.g., Fukue 2000). Since the spectral 
shape of the model could change significantly with the p- value, 
we fitted the XMM1 and XMM2 spectra as well as those from 
Swift and Chandra. The XMM3 spectrum was not considered 
because any evidence for a disk component is only marginal (see 
Servillat et al. 2011). The results are summarized in Table 4. In 
most cases, the p-values are not well constrained. So we cannot 
tell from this model whether or not radial advection plays a role. 

4.2. The Kawaguchi (2003) Disk Model 
4.2.1. Description of the Model 

Kawaguchi (2003) computed disk spectra in a self-consistent 
way, taking into account the effects of electron scattering 
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(opacity and disk Comptonization) and the effects of the 
relativistic correction (i.e., gravitational redshift and transverse 
Doppler shift) on the disk effective temperature in the inner 
part of the accretion disk. Not considering Doppler boosting 
implies that disk solutions are mostly seen face-on ( i = 0 °). 
The model table of the disk spectra , 11 from sub- to super- 
Eddington accretion rates with BH masses of 1— 10 3 Mq, was 
successfully used for analysis of several ULXs (e.g., Foschini 
et al. 2006; Vierdayanti et al. 2006; Okajima et al. 2006; Yoshida 
et al. 2010). We used here a new table including additional 
computation for higher BH masses, extending up to 10 5 Mq. 
The key parameters of the model are the BH mass (M), the 
accretion rate (M), and the viscosity parameter (a). For sub- 
Eddington accretion rates (i.e., in the case of the standard disk), 
the emergent disk spectra are a -insensitive. Near and at super- 
Eddington rates, however, the disk spectra become a-sensitive 
because electron scattering, which is density sensitive (and 
density is a-sensitive), begins to affect the emergent spectra 
quite strongly. The normalization is fixed using the source 
distance ( d — 95 Mpc). To perform the fits, we considered the 
model option 7 that takes into account the effects of advection, 
electron scattering on opacity, Comptonization, and relativistic 
effects (see Appendix A.l for more details about the influence 
of the different model options on the spectral parameters from 
HLX-1 data). It is important to keep in mind that the K03 code 
automatically computes at each radius how much advection is 
present for a given accretion rate. All the models are computed 
assuming a non-rotating BH. 

4.2.2. Results 

To avoid calibration uncertainties between instruments af- 
fecting the fits (see Tsujimoto et al. 201 1), we decided to fit the 
spectra from XMM-Newton, Chandra, and Swift separately. 

We first fitted the Swift-XKT spectra S i p + 2 p , S^^p, Scro, 
and Scri by tying together the BH mass between them. We 
proceeded in the same way with the viscosity parameter a, 
while the accretion rate was left as an independent parameter 
between the different spectra. For the Scro and Sqri spectra, 
we added a power-law component of which the photon index 
was tied together between the two spectra. This is because we 
found consistent values when fitting the two spectra individually. 
When we fixed the a-value to the default K03 model (a = 0.01) 
we did not obtain a good fit because some spectra display strong 
residuals. In order to reproduce the spectra, it was necessary 
to leave the viscosity parameter a free. We then obtained a 
good fit (/ 2 /dof = 108/104). The best-fit parameters are given 
in Table 5. The derived BH mass is M = 1 .8/° 0 2 | x 10 4 Mq. 
Fitting the data from the outburst P3 using the K03 disk model, 
we found M = 1 ,9/q 9 2 x 10 4 M 0 . These estimates are consistent 
and well inside the IMBH mass range. We note that the power- 
law component was not needed to obtain a good fit for the S t 3 
spectrum when leaving the viscosity parameter free to vary. 
Leaving the /V H - value free to investigate the sensitivity of the 
BH mass with respect to the Ah value did not change the 
spectral parameters much and they are still consistent within 
the errors with the values derived with A H = 4 x 10 2 ° cm -2 . 
Indeed, we found A H = 6.5 * 2 8 5 x 10 2 ° cm ~ 2 and M = 
2 . 8*09 x 10 4 M q from the outbursts P0-P2, while we found 

A h = 4. 1 * 2/7 x 1020 cm_2 and M = 1 - 9 -o 2 6 x 1q4 m o from 
the outburst P3. 


1 1 http://heasarc.nasa.gov/xanadu/xspec/models/slimdisk.html 
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Second, we fitted the Chandra spectrum. For the K03 disk 
model, we estimated for the pileup model a grade morphing 
parameter of 0.33. Note that this parameter does not have a 
strong impact on the spectral parameters derived here. The 
derived BH mass is M = 1.9 ± 0.2 x 10 4 Mq. We then left 
the An-value free and we found Ah = 3.7* 25 x 10 2 ° cm -2 and 
a BH mass estimate (M — 1 . 8*04 x 10 4 Mq) consistent with 
those derived from Swift. 

Even if Farrell et al. (2009) showed that the XMM1 spectrum 
is well fitted by a steep power law (T ~ 3.4), suggesting that the 
source was in the steep power-law state as seen in some GBHBs, 
we decided to fit the XMM1 spectrum as well as the XMM2 
one with the K03 disk model. This is because the emergent 
spectrum from the K03 disk model can have a shape strongly 
different from that of a simple multi-color disk spectrum (see 
Appendix A.l). We tied the BH mass between the two spectra, 
but we left a and M to vary independently. First, we keep Ah 
fixed at 4 x 10 2() cnT 2 . For XMM2, it was necessary to add a 
power-law component to obtain a good fit, while the addition 
of a power law for XMM1 did not improve the fit. The derived 
BH mass is then M = 1.4 ± 0.1 x 10 4 M e , which is not 
consistent with the Swift-XKT and Chandra estimates within 
the 90% errors. If we consider the 3a errors, then the estimates 
from the different instruments are all consistent within the errors. 
Then, in order to investigate the sensitivity of the BH mass with 
respect to the Ah - value, we left the Ah parameter free. We 
obtained a value of Ah = 5.5/ 1 ,' l () x 10 2 ° cm ~ 2 and a BH mass 
of M — 1 - 9*0 3 x 10 4 Mq. The photon index of the power-law 
component in both cases is consistent within the 90% errors 
between the XMM-Newton and Swift-XKT spectra. The best-fit 
results are given in Table 5. 

For XMM1, the accretion rate (m ~ 4.4) and a luminosity 
that is very different from that derived using a power law 
(L ~ 1.3 x 10 42 erg cm 2 , a luminosity which was observed 
at the peak of the three outbursts) are commensurate with the 
values found for the XMM2 and Scro spectra (see Table 5). 
This, however, implies a large value of a = 0.13 ± 0.07, 
even if all a-values are consistent within the 3a errors. If 
we force the viscosity parameter to the lowest possible value 
of the K03 model (a = 0.01) keeping the mass fixed at 
M = 1.8 x 10 4 Mq and Ah free, we did not obtain a good 
fit (x 2 /dof = 47.3/24). The addition of a power law improves 
the fit (x 2 /dof = 23.3/22), but in this case the power-law 
component dominates over the disk component. In this case, we 
found similar Y - and A H -values to those given in Farrell et al. 
(2009). 

5. DISCUSSION 
5.1. Constraints on the BH Mass 

To be able to derive a dynamical measurement of the BH 
mass in HLX-1 will be very challenging given the distance of 
the source. So we have to rely on indirect estimates such as 
spectral fitting of the X-ray data. Zampieri & Roberts (2009) 
showed that the BH mass estimate derived using spectral fitting 
can be highly variable depending on the disk model used (see 
their Table 2 ). In the case of HLX-1, Davis et al. (2011) 
used their advanced relativistic disk model to fit spectra with 
different shapes. They found a BH mass within the IMBH range 
(from 3000 Mq to 3 x 10 5 Mq) with extreme and opposite 
assumptions. In this paper, we used the K03 disk model to fit 
X-ray spectra of HLX-1 when the source was in various spectral 
states. We note that this model includes spectra for a wide range 
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Table 5 

Summary of the Spectral Parameters when Fitting the XMM-Newton, Swift, and Chandra Data Using the Kawaguchi (2003) Disk Model 


(1) 

Spectrum 

( 2 ) 

A h 

(3) 

M 

(4) 

M 

(5) 

a 

( 6 ) 

r 

(7) 

Nr 

( 8 ) 
7^ tot 

(9) 

7^ disk 

(10) 

X 2 /dof a 

XMMl 

XMM2 

XMMl 

XMM2 

4 

4 

5 5 +11 
J * -1.0 

c C+l.l 
* 3 * — 1.0 

1.4 ± 0.1 
1.4 ±0.1 

1 0+1-3 
1 ’ — 0.3 
1 9+1-3 
1 >y -0.3 

4.8 ±0.3 

5.9 ±0.4 

4 4+0.3 
*•*-0.5 

5 4+0.5 
D -*-0.7 

q 11 +O.O 8 
u - 11 -0.05 

0 mo +0002 

U.U 1 U _ 0000 

0.13 ±0.07 
0 oio +0 005 

u.uiu _ 0000 

1 a+0.3 
L 6 -0.4 

1 q+0.4 
1 *0-0.5 

1.0 ±0.4 

i q+0.4 
1 ■ — 0.6 

4 2 + ^ -2 
** z -0.6 

5 1+ 0 - 1 
d * a -0.3 

4 9+0-2 
*■"— 0.6 

5 7 +0-1 
d ’ / -0.5 

4 2 +0 - 2 
*- z — 0.6 

4 1+0.1 
*- J - 0.2 

4 9+02 
*■"— 0.6 
4 9+01 

*■"—0.3 

283/215 

283/215 

277/214 

277/214 

Chandra 

Chandra 

4 

3 7 +3 ' 9 
—2.5 

1.9 ±0.2 

1 8 + 1-0 
L °-0.4 

10.2+°^ 

10 - 2 -2 7 7 

0 025 +0 009 
u.uz.j _ 0015 

0 023 +0017 
u.uz.j _ 0013 



1 0 7+0 - 3 
lu - / -0.4 

10 5 +0-5 
lu - J -0.6 

1 0 7+0 - 3 
lu - / -0.4 

10-5 + _°o 5 6 

38/28 

37.6/27 

Slp+2p 
S after! p 

ScRO 

SCR3 

4 

4 

4 

4 

i c+0.2 
0.1 
i c+0.2 
A *°-0.1 
i o+0.2 

A -°-o.i 

i c+0.2 
1 *0-0.1 

lal -2 3 

tO.S/Pifi 
5.5 ±0.5 
3.8 ±0.4 

0 O5 +0 04 
u - UJ -0.02 

0 015 +0 005 
u.uij _ 0005 

0 O 2 O +0009 
u.uzu _ 0006 

0.01 fixed 

1 -7+0.6 

^'-o.s 

i " 7 + 0.6 

^'-o.s 

i 4 + 1.6 
1 ‘*-0.9 
1 a+1.5 
L 6 - 1.0 

1 1 0 +0 - 4 
1 A • — 0.7 

9 9+4-2 

y * y -0.3 

6 2 +0-3 
°- z -0.7 

"> O+0.6 

H-0-o 4 7 

9 9+42 

y - y -0.3 
5 1 + 0-2 

2.3 ±0.5 

108/104 

108/104 

108/104 

108/104 

ScRl 

ScR 2 

4 

4 

1 .8 fixed 
1 .8 fixed 

7.1 ±0.9 
4 q+0.3 

0 0 1 6 +0 028 
u.uio_ 0006 

0 01 9 +0 024 
u.uiy_ 0009 

1 6 +0 - 9 
1 ■ — 1.3 

1 6 +09 
a *°-1.3 

1 3 +2 - 5 
1 * J -1.0 

1 3 +2 - 5 
1 ** 5 -l.l 

8-0-o 5 9 

5-5 + _°o 4 6 

6 9 +0 - 4 

4 4+0.4 
*•*-0.5 

46/50 

46/50 

5.1 

5.2 

S, 3 

4 

4 

4 

1.8 fixed 
1.8 fixed 
1.8 fixed 

10.92 0 ! 3 

9.1 ± 1.0 

fi 7+ 1 - 4 
0/ — 1.2 

0 01 +011 
u - ul -o.oo 

0 016 +0 005 
u.uio_ 0006 

n m q+0.023 
U.Uio_o 007 



10.8 ± 1.1 
8.8 ±0.6 
6.7 ±0.9 

10.8 ± 1.1 
8.8 ±0.6 
6.7 ± 0.9 

37.7/41 

37.7/41 

37.7/41 


Notes. Columns: (1) spectrum name; (2) absorption column in units of 10 20 cm~ 2 ; (3) BH mass in units of 10 4 (4) accretion rate in units of l-rdd/c 2 with 

L p,dd and c, the Eddington limit and the speed of light, respectively; (5) viscosity parameter; (6) photon index of the power-law component; (7) normalization 
of the power law in units of 1 0 2 photons keV 1 cm~ 2 s _1 ; (8) unabsorbed 0.3-10 keV total luminosity in units of 10 41 erg s -1 ; (9) unabsorbed 0.3-10 keV 
disk luminosity in units of 10 41 erg s -1 ; (10) / 2 -value and number of degrees of freedom. 

a The Swift-XRT. XMM-Newton EPIC-pn, and Chandra /ACIS data were fitted separately. For a given instrument, we fitted together all the available spectra or 
a subset of them. 


of accretion rates from sub- to super-Eddington rates. In all cases 
for ULXs, the fits favored super- or near-Eddington accretion 
onto a stellar mass BH. However, the fits to HLX-1 spectra do 
favor an IMBH solution with a BH mass of M = 1 .8 +l (| ^ x 
10 4 Mq. This is an interesting result because we get three 
independent and nonetheless consistent mass estimates from 
three different instruments (Swift-X RT, XMM-Newton EPIC-pn, 
and Chandra/ ACIS) and over different spectral shapes. This 
estimate is consistent with the observational lower and upper 
limits on the BH mass (Farrell et al. 2009; Servillat et al. 201 1). 
So, using the K03 model which assumes a non-spinning BH 
and a face-on accretion disk, we find that an IMBH solution is 
favored with an estimated mass of ~ 1.8 x 10 4 M e . The BH 
accretes at the Eddington limit ( m ~ 10) and radiates close 
to the Eddington luminosity at the outburst peak (Ledd = 
2-3-06 x 10 42 er 8 s_1 )- This corresponds to an Eddington 
fraction of /e dd = 1.1*0® considering a bolometric luminosity 
of 2.5 x 10 42 erg s -1 on average (see Figure 5). 

5.2. Comparison with the Results from Davis et al. (2011) 

Using their BHSPEC model, Davis et al. (2011) studied 
the dependency of the BH mass in HLX-1 on the inclination 
and the BH spin. We can only compare our results with their 
work for a* = 0 and i = 0°. In that configuration, they 
found a consistent mass estimate for XMM2 and Chandra with 
log M ~ 3.8 ± 0.1 (S. W. Davis 2011, private communication). 
However, no solution could be found in this case for the Swift 
spectra, because the /ndd = Lx /L em ratio was pegged to 
the maximum permitted value in their BHSPEC model (i.e., 
/kdd = 1; S. W. Davis 2011, private communication). From 
Figures 2, 5, and 6 in Davis et al. (201 1), the best-fit contours for 
fixed inclination in the a* versus log M plots overlap well with 



Accretion rate (L Edd /c 3 ) 


Figure 5. Linear correlation between the 0.01 and 20 keV unabsorbed luminosity 
derived using the Kawaguchi (2003) slim disk model and the accretion rate. The 
luminosity was computed using the following relation: Lbol = 4jzdj Fbol with 
Fbol, the bolometric unabsorbed flux and di = 95 Mpc. 

(A color version of this figure is available in the online journal.) 

the simulations for the range of allowed parameters. Thus, they 
found a consistent mass estimate between the three instruments 
for i = 0° and a* = 0.7 with log M ~ 4.0 for Swift and Chandra 
and log M ~ 3.9 for XMM2. The mass values derived by Davis 
et al. (2011) from the XMM2 and Chandra spectra for a* = 0 
and i = 0° are smaller than the one derived using the K03 disk 
model. We note that the Eddington ratios (/Edd) they obtained 
are large (and close to 1). From our fitting results using the 
K03 model, we found similar high /'em values near the outburst 
peak. We believe that for such high /Edd values the effects of 
advection on the emergent spectra should be considered. 
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Evaluating the level of systematic errors of the K03 model 
by comparison with alternative models by exploring the whole 
parameter space is beyond the scope of the present paper and tr 

deserves a dedicated paper. Nevertheless, we estimated the level \ 

of systematic errors of the K03 model by comparison with the jf 

KERRBB (Li et al. 2005) and BHSPEC models assuming a ~ ' 

non-spinning BH ( a * = 0) and a face-on disk (i = 0°) in — 

the limit of low mass accretion rate (i.e., a few percent of the s- 

Eddington limit). In that limit, the accretion disk is expected fl 

to be very close to the relativistic model of Novikov & Thorne 
(1973). We found that the level of systematic errors in the BH 
mass between the K03 and KERRBB models is less than 45% 
when only considering the relativistic effects in both models (see 
Appendix A. 2). When comparing the results between the K03 
and BHSPEC models taking into account the relativistic effects 
and the effects due to Comptonization and electron scattering, 
we found that the level of systematic errors in the BH mass is 
less than 24% (see Appendix A. 2). In all cases, the values of the 
accretion rate derived using the K03, KERRBB, and BHSPEC 
models agree within the errors. So, providing we have a similar 
level of systematic errors over the whole parameter space, we 
are confident that our BH mass and accretion rate estimates are 
reliable. 

5.3. The Disk Structure 

Figure 5 shows a linear correlation between the 0.01-20 keV 
bolometric disk luminosity (L) and the accretion rate so that 
L = fph L|-dd with radiation efficiency r] = 0.1 1 ±0.03. Watarai 
et al. (2000) found a similar linear correlation for m < 20. This 
implies that advection does not play a significant role, and that 
the slim disk solutions are not needed. At all times during 
the outbursts, the disk aspect ratio H/R is less than 0.2, from 
Figure 5 in K03. Note that H/R is mostly independent of the 
BH mass and the viscosity parameter. 

From Figure 1 in K03, the disk appears to be radiation 
pressure dominated in its inner regions for all the m-values 
found. Such disks might be viscously and thermally unstable, 
and they might give rise to outbursts (Taam & Lin 1984; Honma 
et al. 1991; Lasota & Pelat 1991; Xue et al. 2011). However, 

Hirose et al. (2009) showed via MHD simulations that such 
disks could be thermally stable (but see Xue et al. 201 1). Most 
observations of Galactic X-ray binaries in outburst radiating 
above 10% of the Eddington luminosity for a stellar mass BH 
do not show any evidence for such an instability except perhaps 
for GRS 1915+105 (Belloni et al. 1997; Xue et al. 2011). 

5.4. The Evolution of the Accretion Rate 

Figure 6 shows the evolution of the accretion rate through the 
different outbursts. At the peak of the three outbursts and for a 
few weeks (~3 and ~2-3 weeks for the outbursts P2 and P3, 
respectively), the BH accretes at the Eddington limit (m ~ 10). 

This “plateau” is seen in the X-ray light curve of the outbursts 
P2 and P3 (see the right panel in Figure 1). Recently, Lasota 
et al. (201 1) proposed that the HLX-1 outbursts may be due to 
enhanced mass-transfer rate onto a pre-existing accretion disk 
when an asymptotic giant branch star orbiting along an eccentric 
orbit with a period of ~ 1 year passing at periastron is tidally 
stripped. Such a mechanism is known to produce “plateaus” 
in light curves (Bath & Pringle 1981; Esin et al. 2000). The 
latter authors considered an unstable disk that became stable 
after an enhanced mass-transfer rate phase. In the case of 
HLX-1, the situation is different since the disk is stable. 



54800 55000 55200 55400 55600 55800 
Time (MJD) 

Figure 6. Evolution of the accretion rate over time. 

(A color version of this figure is available in the online journal.) 

Therefore, a quasi-constant enhanced mass-transfer rate during 
a given time interval would naturally produce a “plateau” in 
the light curve. After the “plateau,” the accretion rate drops 
exponentially with a decay time of 182 ± 30 days and 92 ± 15 
days for the outbursts PI and P2, respectively. The difference in 
the decay time is because there was a small re-flare event at the 
end of the outburst PI (see Figure 1). From the evolution of the 
accretion rate over time, we calculated the mass accreted during 
the outbursts PI and P2 assuming that the accretion rate is the 
same throughout the disk. If this is not the case, then the values 
derived below are lower limits. We found an accreted mass of 
~1 .2 x 1 0 4 M q and ~8.2x 10 -5 M 0 for the outbursts PI and 
P2, respectively. 

At the end of the outbursts PI and P2, i.e., when the source 
undergoes a transition from the high/soft state to the low/hard 
state, the count rate drops very quickly over a timescale less 
than a month (see Figure 1). As shown by Servillat et al. 
(2011), any disk component in the spectrum of HLX-1 in the 
low/hard state is marginal. Given the XMM3 luminosity of 
~2 x 10 40 erg s -1 derived from the power-law component and 
assuming a radiation efficiency in the low/hard state of 0.11 
as during the outbursts, we could compute an upper limit on 
the accretion rate in the low/hard state of m < 0.09, i.e., 
M <3.6 x 1 0 6 M !:) yr 1 . This is nicely consistent with the 
results obtained by Esin et al. (1997, 1998) using the advection- 
dominated accretion flow model. According to that model, this 
must be connected with a varying (receding) inner disk radius 
(e.g., Dubus et al. 2001). Higher statistical quality of the data 
when the source transits to the low/hard state will be needed to 
further investigate that possibility. 

5.5. The L-T Relation 

In the luminosity-temperature space, we found a correlation 
L oc T~ 2A . Fukue (2000) and Kawaguchi (2003) showed that 
near the Eddington limit (hi = 10) and above a departure 
from the L-T 4 relation is expected, due to the change of the 
radial profile of the disk temperature and the m -sensitive disk 
color temperature via the effects of electron scattering. However, 
below this limit the disk is still expected to follow the LcxT 4 
relation found in the standard disk model. Apart from the peak 
of the outbursts when the accretion rate is at the Eddington limit, 
the accretion rate is below this limit. Therefore, we should expect 
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to have a LocT 4 correlation. The behavior near the peak is 
probably responsible for the flatter correlation. To check that, we 
fitted our L-kT points for the XMM2, .S'cr i , 5 cr 2, ■'>CR3, and S , 3 
observations (for which m < 10 — see Table 5) using the relation 
for a standard disk. We found: log(L) = 4 log T + 44.95 ± 0.06 
(la error). From Equation (1) in Lasota et al. (2011) giving 
L as a function of the color disk temperature at the inner disk 
radius and the BH mass, we derived the following relation: 
log(L) ~ 41ogT +45.1 using ;; ~ 0.11, M = 1.8 x 10 4 Mq, 
and /?isco = 3 Rs (for a non-rotating BH). This demonstrates 
that the disk in HLX- 1 follows the L oc T A relation. 

6. CONCLUSION 

We investigated in detail the X-ray spectral properties of the 
best IMBH candidate, HLX-1, using multi-epoch data collected 
by three different X-ray instruments (Swift-X RT, XMM-Newton 
EPIC-pn, and Chandra). To do so, we used the K03 disk model 
in order to constrain the BH mass, the accretion rate, and the 
disk structure in a non-biased way. Indeed, this disk model 
has the advantage of taking into account the effects of radial 
advection and covering a wide range of accretion rates (from 
sub- to super-Eddington) assuming a non-rotating BH and a 
face-on disk inclination (i = 0°). We found that our multi- 
epoch data are consistent with sub-Eddington accretion onto a 
nearly 2 x 10 4 M e BH. At the peak, the X-ray luminosity is 
near the Eddington luminosity. The derived radiation efficiency 
is r) ~ 11%. The disk solution we found for HLX-1 is different 
from those derived for other ULXs using the same disk model 
(e.g., Vierdayanti et al. 2006; Okajima et al. 2006; Yoshida 
et al. 2010). Indeed, for other ULXs the spectral fitting favors 
super- or near-Eddington accretion onto a stellar mass BH. Here, 
the disk in HLX-1 likely undergoes sub- and near-Eddington 
accretion onto an IMBH and appears to stay relatively thin. The 
inner regions of the accretion disk are dominated by radiation 
pressure. We also showed that the source follows the L oc T 4 
relation globally. At the outburst peak and for a further few 
days, the BH accretes at near the Eddington limit (m ~ 10, i.e., 
~4 x 1 0 4 M q yr ' ). The occurrence of this “plateau” at the 
outburst peaks could provide some evidence that the outburst 
mechanism is driven by enhanced mass-transfer rate onto a pre- 
existing accretion disk as proposed by Lasota et al. (201 1). After 
the “plateau,” the accretion rate decreases exponentially until the 
end of the outburst. However, it is not clear from our data what 
is the geometry of the accretion flow when the source transits 
to the low/hard state, even if the presence of an optically thin 
advection-dominated flow would avoid a dramatic variation in 
the accretion rate by more than two orders of magnitude from the 
peak to the low/hard state. Finally, the recent outburst starting 
in 2011 August gives more credit to the ~ 1 year recurrence 
timescale proposed by Lasota et al. (2011). Our dedicated 
multi-wavelength (radio, optical, and X-ray) observations of 
this outburst will shed further light on the nature of this unique 
source as well as on the outburst mechanism. 
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Science and Technology Funding Council. S.A.F. is the recipi- 
ent of an Australian Research Council Postdoctoral Fellowship, 
funded by grant DPI 10102889. V.B. acknowledge support from 
the UK STFC research council. J.P.L. acknowledges support 
from the French Space Agency CNES. 

APPENDIX 

A. /. Impact of the K03 Disk Model Computing Options 
on the Shape of the Emergent Spectra, the BH Mass, 
and the Accretion Rate 

The K03 model makes use of different options to compute the 
emergent spectra starting from a local BB and then including 
different physical effects (advection, electron opacity, relativis- 
tic effects, and Comptonization). First, we investigate the impact 
of these different options on the shape of the emergent spectra 
as well as on the M and M parameters. We take the example 
of the S\ p + 2 p spectrum as an illustration, but what is described 
below is also valid for the other spectra. 

Figure 7 shows the modification of the S\ P+ 2 P spectrum 
when selecting different model options: (1) standard disk and 
local BB; (2) slim disk and local BB; (5) slim disk including 
relativistic effects on the BB; (6) slim disk including the 
effects of electron scattering on opacity and relativistic effects; 
(7) slim disk including the effects of electron scattering on 
opacity, Comptonization, and relativistic effects. Here, the term 
“slim disk” implies that the K03 code automatically computes 
(at each radius) how much advection is present for a given 
accretion rate. 

First, we just changed the option parameters without running 
a new fit for option 1 in order to see the effects of the 
different model components. When using option 1, we got 
for the S\ p+ 2 p spectrum a BH mass around 1.8 x 10 3 Mq and 
m ~ 32 assuming a viscosity parameter of a = 0.01. Such high 
m -values mean super- Eddington accretion. However, in this case 
radial advection is expected to start playing a significant role. 
The large derived M-value also implies that electron opacity 
and Comptonization should play an important role (see Figure 2 
in Kawaguchi 2003, for instance). This results in more spectral 
boosting toward high energies. In order to obtain the global 
shape of the data when including all the above effects in a 
slim disk (option 7), it is clear that the effects of both electron 
scattering on opacity and Comptonization have to be reduced. 


10 


The Astrophysical Journal, 752:34 (12pp), 2012 June 10 


Table 6 

Evolution of the Parameters of the Kawaguchi (2003) Slim Disk Model for 
Different Model Options 


Option 

(i) 

M 

(2) 

m 

(3) 

a 

(4) 

X 2 /dof 

(5) 

1 

2 

5 

6 
7 

i q+0.4 
0.2 
n r+0.3 
D.D—0 2 

r-) 0 + 1. 0 
z " 3 — 0.6 
90+4 

” — 14 
18 ± 2 

32.27 1 5 1 6 

38.02 1 ! 4 

200 +477 

zuu -101 

lO.tP.o.9 

10-1-23 

0.01 fixed 

0.82$$ 

>0.26 

0 iH+O-Oi 
u - 1 ^-0.04 

0 05 +0 06 
u,u - 0.02 

40.4/36 

40.1/35 

51.59/35 

67.21/35 

43.4/35 


Notes. The fits were performed using the Si p + 2 P spectrum with Nu = 
4 x 10 20 cm~ 2 . Columns: (1) computation option from the Kawaguchi (2003) 
slim disk model (see Section A.l); (2) BH mass in units of 10 3 M Q ; (3) accretion 
rate in units of L Edd/c 2 ; (4) viscosity parameter; (5) value of the y 1 and number 
of degree of freedom. 

To do so, the accretion rate (in units of the Eddington rate) has 
to decrease, which in turn results in an increase of the BH mass. 

To check this, we fitted the S\ P+ 2 P spectrum for each chosen 
option leaving the model parameters free. The evolution of the 
model parameters are summarized in Table 6. From option 1 
to option 2, there is a relatively small variation of the M and 
m parameters (M ~ 3.4 x 10 3 M 0 and m ~ 38). From option 
2 to option 5 (including the relativistic effects), we obtained a 
significant impact on the accretion rate (m ~ 200) and a slight 
decrease of the BH mass (M ~ 2.3 x 10 3 Af G ). Given the large 
accretion rate (m ~ 38) derived in option 2, the accretion disk 
extends below the last stable orbit (e.g., Watarai et al. 2000). The 
relativistic effects (gravitational redshift and transverse Doppler 
redshift) when applying option 5 strongly suppress the emission 
coming from these regions, which in turn results in the model 
underestimating the data. To have a good match between the 
data and the model, the accretion rate has to be increased. From 
option 5 to option 6, the spectral parameters significantly change 
with an increase in the BH mass by one order of magnitude and 
a decrease of the accretion rate by more than one order of 
magnitude. From option 6 to option 7, the parameters do not 
change significantly (M ~ 1.8 x 10 4 M 0 and m ~ 10). 

From option 5, the accretion flow is in the super-Eddington 
regime. The effects of Comptonization and electron scattering 
(option 6) are likely to be significant given the large m -value 
(see Figure 2 in K03). This would in turn result in a stronger 
hardening of the model than what is seen in Figure 7. However, 
in order to obtain a good fit given the shape of the S i p + 2 p 
spectrum, the effects of both Comptonization and electron 
scattering have to be strongly minimized. This is achieved by 
significantly decreasing the accretion rate (the accretion flow is 
then in the sub-Eddington regime) and therefore by significantly 
increasing the BH mass. If we force the source distance from 
95 Mpc (measured distance) to an arbitrary lower value (e.g., 
cl = 3.5 Mpc), such a dramatic change in the BH mass and 
the accretion rate is not seen because the accretion flow always 
stays in the sub-Eddington regime (see Table 7). 

A. 2. Comparison of the K03, KERRBB, and BHSPEC 
Results in the Limit of the Novikov & Thome (1973) 
Relativistic Disk Model 

We compare the results from the K03 disk model with those 
obtained with KERRBB and BHSPEC assuming a non-spinning 
BH (a* = 0) and a face-on disk (i = 0°) in the limit of low 
mass accretion rate (i.e., less than 10% of the Eddington limit). 


Godet et al. 

Table 7 

Evolution of the Parameters of the Kawaguchi (2003) Slim Disk Model for 
Different Model Options Forcing the Source Distance to 3.5 Mpc Instead of 
the Measured Distance of 95 Mpc 


Option 

(i) 

M 

(2) 

m 

(3) 

O' 

(4) 

X 2 /dof 

(5) 

1 

2 

5 

6 
7 

70+15 

/v -12 

642\ 5 o 

587‘s 4 

52+ 3 

122 +13 

iZZ -10 

1 07+0-05 
A - u/ — 0.04 

1 7+0- 2 
A * —0.2 
0 q+0.4 
z ‘ — 0.5 

0 O+ 0. 1 

— 0.1 

1 4+0.1 
A *^-0.1 

0.01 fixed 
0.01 fixed 
0.01 fixed 
0.01 fixed 
0.01 fixed 

40.2/36 

40.3/36 

40.8/36 

40.4/36 

42/36 


Notes. Columns: (1) computation option from the Kawaguchi (2003) slim disk 
model (see Section A.l); (2) BH mass in units of M 0 ; (3) accretion rate in units 
of f-Edd/c 2 ; (4) viscosity parameter; (5) value of the y 2 and number of degree 
of freedom. 

In that limit, the accretion is expected to follow the predictions 
of the relativistic disk model of Novikov & Thorne (1973). In 
the K03 model, this limit corresponds to m f 1 .('iLppy/c 2 . For 
the comparisons with KERRBB, we limit ourselves by setting 
the hardening factor to 1 in the model (i.e., we did not take the 
effects of electron scattering and Comptonization into account). 
We also did not consider the effects of limb darkening and self- 
irradiation since they are not included in the K03 disk model. 
This configuration corresponds to option 5 in the K03 disk model 
(only including the relativistic effects). For the comparisons 
with BHSPEC, we used the option 7 (including the effects for 
Comptonization and electron scattering) in the K03 disk model. 
We used the table of BHSPEC models bhspec_mass_0.01.fits 12 
considering BH mass up to 300 M 0 . 

To perform the comparisons, we proceeded in both cases as 
follows. 

1. We fitted the S\ P+ 2 P spectrum (i.e., when the source lu- 
minosity is maximum) using the WABS*KERRBB and 
WABS*BHSPEC model with a* = 0, i = 0°, and 
Nu = 4 x 10 20 cm~ 2 for different and arbitrary values 
of the source distance D so that the m-values derived were 
below 10% of the Eddington limit. 

2. From the best fits, we simulated a model spectrum with an 
exposure time of 38 ks using the Swift-X RT response files 
and the fakeit command in Xspec. 

3. We fitted the faked spectra with the WABS*K03 model 
with a fixed to 0.01, the same source distance D as used 
for the KERRBB and BHSPEC fits and the options 5 and 7 
when comparing the results with KERRBB and BHSPEC, 
respectively. 

Tables 8 and 9 summarize the results of the comparisons 
with KERRBB and BHSPEC, respectively. From Table 8, we 
show that there is a systematic discrepancy (less than 45% in 
the range considered) in the BH mass between the two models, 
while the accretion rate values appear to be consistent within the 
errors. From Table 9, we obtain a good agreement between the 
BHSPEC and K03 models within the 3 ct errors except for the 
BH mass estimate for a distance of 2.9 Mpc. Based on that latter 
case, the “systematic errors” in the BH mass are less than 24%. 
We also show in Tables 8 and 9 comparisons with hi -values 
larger than 10% of the Eddington limit. Again, we emphasize 
the good agreement between the K03-KERRBB models and 
K03-BHSPEC models. So, we conclude that the K03 disk model 


12 http://www.cita.utoronto.ca/~swd/xspec.html 
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Table 8 

Summary of the Results Obtained for both the KERRBB and K03 Models 


Distance 

Parameter 

KERRBB a 

K03 b 

m c 

(Mpc) 




(^Edd/ 

10. 

Mass (M 0 ) 

145 

19U 4 2 2 9 d 


10. 

m (10 19 g s _1 ) 

21.6 

18-8 ^ 

( 6 - 82 . 0 . 9 ) 

5. 

Mass (Mq) 

73 

105* g 


5. 

m (10 19 g s _1 ) 

5.4 

r O+0.4 
— 0.4 

(3-SS.‘i) 

2.9 

Mass (Mq) 

42 

522 9 7 


2.9 

m (10 19 g s _1 ) 

1.8 

1 C ) +0A 
A ' — 0.3 

(2.02°o 3 2 ) 

2. 

Mass (Mq) 

29 

39* 4 


2. 

m (10 19 g s -1 ) 

0.85 

0.75 ± 0.08 d 

(1.322 0 o°o 7 7) 

1.5 

Mass (Mq) 

22 

32*j 


1.5 

m (10 19 g s -1 ) 

0.49 

0 50 +0 04 
u.JU—008 

(1.0B + _°o°o 4 3) 


Notes. The errors are given at the 90% confidence level. 

a The effects of Comptonization, electron scattering, limb darkening, and self- 
irradiation were not taking into account. 

b We used the option 5 in the K03 model that only includes the relativistic 
effects. 

c Accretion rate value derived from the K03 disk model. 
d Consistent at the 3er level. 


Table 9 

Summary of the Results Obtained for both the BHSPEC and K03 Models 


Distance 

(Mpc) 

Parameter 

BHSPEC 

K03 a 

m° 

^Edd/ ^ 

7. 

7. 

Mass (Mq) 

l = Lx /L eu 

234 

0.146 

2372 7 2 3 i 

0.2002°o°o 6 4 C 

(3.202°o 0 o 6 4 ) 

5 . 

5 . 

Mass (Mq) 
1 

158 

0.110 

172 +31 

1 23 

nil £+ 0.024 
0 . 116_0 020 

(1.92°o 2 2 ) 

2.9 

2.9 

Mass (Mq) 
1 

85 

0.069 

1052g d 

0 07 3+0.006 
U.U / j_q 006 

d-l® 

2. 

2. 

Mass (Mq) 
1 

54 

0.051 

5521 

0 064 +0 007 
U.UO‘+_o 006 

a.02*_° 0 °02 e ) 


Notes. The errors are given at the 90% confidence level. 
a We used the option 7 in the K03 model that includes the relativistic effects as 
well as the effects of Comptonization and electron scattering. 
b Accretion rate value derived from the K03 disk model. 
c Consistent at the 3 a level. 
d Marginally consistent at 3er . 

e Lower error bar pegged to the lowest permitted value of the model (m = 1). 

provides reliable estimates of accretion disk parameters (mass 
and accretion rate). 
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